Enhancing multiphoton rates with quantum memories 



J. Nunn^'QN. K. Langford, 2 W. S. Kolthammer, 1 T. F. M. Champion, 1 
M. R. Sprague, 1 P. S. Michelberger, 1 X.-M. Jin, 1 ' 3 D. G. England, 1 and I. A. Walmsley 1 

1 Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom 
Department of Physics, Royal Holloway, University of London, 
Egham Hill, Egham TW20 OEX, United Kingdom 
1 Centre for Quantum Technologies, National University of Singapore, 117543, Singapore 

(Dated: August 9, 2012) 

Single photons are a vital resource for optical quantum information processing. Efficient and 
deterministic single photon sources do not yet exist, however. To date, experimental demonstrations 
of quantum processing primitives have been implemented using non-deterministic sources combined 
with heralding and/or postselection. Unfortunately, even for eight photons, the data rates are 
already so low as to make most experiments impracticable. It is well known that quantum memories, 
capable of storing photons until they are needed, are a potential solution to this 'scaling catastrophe'. 
Here, we analyze in detail the benefits of quantum memories for producing multiphoton states, 
showing how the production rates can be enhanced by many orders of magnitude. We identify 
the quantity r\B as the most important figure of merit in this connection, where r\ and B axe the 
efficiency and time-bandwidth product of the memories, respectively. 

PACS numbers: 42.50.Ex, 42.50.Ct, 42.50.-p 



After two decades of rapid advances, quantum op- 
tics experiments are becoming increasingly challeng- 
ing. As the interests of the community shift to higher- 
dimensional entanglement [T] [2] and information process- 
ing tasks beyond mere proof-of-principle [3l H], the de- 
mand for large numbers of simultaneous single photons is 
outstripping the capabilities of parametric sources [51 [5] . 
These sources, which so far have been the workhorse of 
the quantum optics lab, produce photons in pairs, but 
they also produce multiple unwanted photon-pairs with 
a probability that scales with the single-pair generation 
rate, which must therefore be kept low, so that most of- 
ten no photons are emitted. The current record for pho- 
tonic resources is an eight-photon experiment involving 
four parametric sources, in which statistics were accu- 
mulated over 40 hours [7]. Besides the exorbitant time 
required to run experiments of this kind, the quality of 
the multiphoton states produced is limited by multi-pair 
emission. One solution is to operate the sources in a 
heralded fashion, in which one of each photon pair is de- 
tected, indicating the presence of the other photon. This 
halves the number of usable photons but it removes the 
non-emission events, which makes noise events a much 
smaller fraction of the finally-detected photons. The use 
of heralded parametric sources is extremely convenient, 
since they operate at room temperature and at extremely 
high repetition rates. Nonetheless the low success proba- 
bility means they cannot be efficiently combined to pro- 
duce multiphoton states. 

The scalability problem is easily understood. Suppose 
that a single photon is heralded with probability q -C 1. 
The probability of producing N single photons simultane- 
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ously using N sources is then simply q N , which becomes 
exponentially small as N increases, thus rendering com- 
plex experiments impossible. 

One method to mitigate this problem is by multiplex- 
ing many sources via active switching [8HT3"]. but large 
overheads — i.e. many identical sources — are required 
to achieve efficient operation [2]. It is well-known that 
temporal multiplexing using quantum memories offers an 
alternative solution [T5HT?] . 

To see how quantum memories can increase the rate 
of iV-fold coincidences, consider the array of N sources 
coupled to N — 1 memories shown in Fig. [T] We suppose 
that each source produces photons in pairs by means of 
a parametric scattering process such as downconversion 
[6] or spontaneous four- wave mixing [5] , with one of each 
pair directed to a herald detector. With no memories, all 
TV heralds must fire simultaneously to produce an A^-fold 
coincidence. However with memories, heralded photons 
can be stored whenever they are produced. Once all N— 1 
memories are charged with a photon, one only has to wait 
for the final source to produce a photon, and then all the 
memories can be read out and one has, again, an A^-fold 
coincidence. This protocol is probably not optimal, but 
it is amenable to a straightforward analysis that captures 
the scaling enhancement: by lifting the requirement for 
simultaneous emission, the memories greatly enhance the 
coincidence probability. Our purpose in this paper is to 
quantify the gain in coincidence rate afforded by using 
quantum memories to synchronize photon sources in this 
way. The time-bandwidth product B = St proves to be 
critical in this context, where S is the acceptance band- 
width of the memories, and r is their coherence lifetime 
[32] . If postselection on the final detection of N photons 
is used, even relatively inefficient memories can dramat- 
ically enhance the multiphoton rate. 

We first fix the protocol by assuming that we always at- 
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each memory can be tracked using a transfer matrix: 



FIG. 1: An array of N heralded parametric sources synchro- 
nized by N— 1 quantum memories. The sources are repeatedly 
pumped, and each photon emitted is stored until all memories 
are charged. Then emission of a photon by the final source 
triggers retrieval from the memories, in order to generate an 
TV-fold coincidence. 



tempt to store a photon, if a herald detector fires, regard- 
less of whether or not the memory concerned is already 
charged. This ensures that we always use the most re- 
cently emitted photons, which mitigates photon loss due 
to decohercncc in the memories. To avoid 'clashes' (in 
fact, interference [HI [19]) between incident and stored 
photons, we clean each memory before storage is at- 
tempted (e.g. by readout of the memory, or optical 
pumping) so that we are always attempting to charge 
an empty memory. This allows to use a classical model, 
in which individual photons are treated as particles that 
are probabilistically emitted, stored and retrieved. Fi- 
nally, we adopt the policy that if we are ready to read 
out the memories — that is, if all memories have been 
charged and a photon is emitted from the N th source — 
and at the same time one or more of the other sources 
emits a photon, we bypass the relevant memories and use 
these 'serendipitous' photons, rather than attempting to 
read out the memories. 

The photon sources are pumped at a rate 1Z ~ <5, lim- 
ited by the minimum pulse duration that can be stored 
by the memories. The average waiting time 1/TZc be- 
tween TV-photon events can then be computed if we can 
find an expression for the TV-fold coincidence probability 
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which is the probability that one photon is 
obtained from each of the N sources. The leading factor 
of q describes the probability that the N th source emits a 
photon, and we have defined p S ync = Q+QV^P &s the prob- 
ability that any one of the N — 1 sources equipped with a 
memory provides a photon on demand, either directly, or 
through successful retrieval of a stored photon. Here P is 
the steady-state probability that any memory is charged 
with a stored photon, ?y r is the retrieval efficiency, and the 
overbar notation denotes the probabilistic complement, 
X = 1 — X. The problem of computing the waiting 
time then reduces to that of finding P. To proceed we 
assume that the decoherence processes in the memories 
are Markovian (i.e. exponential), since then the stochas- 
tic evolution of the charge-state x^ = [P (m) , P< m )]T of 
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with P(™) the probability that the memory is charged at 
the m th time step, r the probability that an empty mem- 
ory becomes charged over the course of one time step, 
and s the probability that a charged memory becomes 
empty. The steady state probabilities are given by the 
eigenvector x s of T with eigenvalue 1, x s = [s, r] T / (r+s), 
so that we have P = r/(r + s). The probability that an 
empty memory becomes charged is the probability that a 
heralded photon is emitted and that it is stored, provided 
that the rest of the set-up is not primed for readout, so 
we have r — qr) s R, where R is the probability that the 
system is ready to be read out (the evaluation of R is de- 
scribed in the Appendix), and r] s is the storage efficiency. 
The loss probability that a charged memory is emptied 
is more complicated. There are four processes involved. 
First, decoherence in the memory during standby, sec- 
ond, readout of the memory when we attempt to generate 
a coincidence, third, the loss of a stored photon during 
standby when a new photon comes along and we attempt 
to replace the stored photon but fail, and finally deco- 
herence in the memory during the readout stage, when a 
photon is heralded and the memory is bypassed, leaving 
the memory charged and vulnerable to decay. Denoting 
the decoherence probability by b, the total loss probabil- 
ity works out to be s = q[bR + R] + q[rj s R + Rb]. Gen- 
erally the time-bandwidth product will be much larger 
than one, so that 6=1 — e~ x l B ss l/£>, and we finally 
obtain 



C = q N il + 



RqrjB 



l + (B-l)[R(q-q) + q] 



N-l 



(2) 



This is the main analytic result of this paper. In the limit 
of small photon generation rates such that {RB,qB} <C 
1, we have c ~ q N (r]B) N ^ 1 , which supports the intu- 
ition that each memory effectively boosts the photon 
generation probability by B, moderated by its efficiency 
rj = ry s ?7 r . In this regime the gain in the multiphoton rate 
is therefore exponential in the quantity rjB, which high- 
lights the importance of the time-bandwidth product for 
synchronization applications. As B is increased so that 
qB 3> 1, the rate eventually saturates and becomes inde- 
pendent of B, limited finally by r\. 

To make a fair comparison with the unsynchronized 
case, we now consider the effect of higher photon num- 
ber components on the quality of the states produced. 
Typically, parametric sources generate photon pairs ac- 
cording to a thermal distribution, where p SOU rce( rl ) = PP n 
is the probability of emitting n photon pairs, and p 
is a small real number. We also assume non-photon- 
number-resolving heralding detectors, such as APDs, so 
that the conditional probability that n photons are sent 
towards a memory, given a herald click, is Ph(n) — 
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[1—h ]j3 S ourcc (n) I q where h is the efficiency of the herald- 
ing detector, and as before q = hp/[l — ph] is the proba- 
bility of a herald click. For simplicity we assume that de- 
tector dark counts are negligible. The charge state 
of the memory is now a vector of probabilities that the 
memory contains n photons, with n = 0, 1, 2..., which we 
truncate for numerical convenience. The transition prob- 
ability that the number of excitations stored in a memory 
changes from k to j over the course of any time step is 
given by the transfer matrix element 



r jk = e jk b k - j b 3 



(Rq + Rq) + qR6 j0 + qR Ps (j), (3) 



where the three terms represent decoherence, readout, 
and storage, respectively. Here djk — 1 for k > j and 
zero otherwise, and Sjo is a Kronecker delta. We have 
also defined p s (n) as the probability that n photons are 
stored in the memory when read-in is attempted after 
a herald, p s (ri) = J2T=n Ph(k)VsVs Repeated ap- 
plication of T to an arbitrary initial charge state con- 
verges to the steady state x s , and the probability that 
n photons are retrieved from the memory is then given 
by Pr(n) = EZn^VrVr^O- We can then write 
c = gp h (l)p sync (l) Ar_1 , where p sync {n) = qp h (n) + qp T (n) . 
This result for c represents only a minor correction to 
Eq. |2]) , but the treatment of multi-pair emissions is im- 
portant for the fidelity calculation below. 

In many photonic networks, successful operations can 
be postselected on the final detection of at least N pho- 
tons. In this case the fidelity of the postselected states is 
the fraction of these which contain 1 photon per mode. 
We normalise this to the number of modes by taking the 
N th root, 
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Here p>n is the probability that the state obtained from 
the memories/sources comprises N photons or more, and 
we have re- written this in terms of p<jv, the probability 
that fewer than N photons in total are emitted, given by 
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Here Sfc(Z) is the I th element of a vector s& containing 
N — 1 real, non-negative integers whose sum is equal to 
k. The summation V] runs over all such vectors. 

If postselection is not used, we consider the fraction of 
readouts that our desired N photon state, with one pho- 
ton per mode, is produced, to obtain T = [c/RY] 1 ^ , 
where RY is the probability that we believe we have pro- 
duced an N photon state (see Appendix). 

Without memories, the TV-fold coincidence rate is 
c no mom = [qPh{l)] N and the unpostselected fidelity is 
simply given by J-" no mem = Ph(l). The postselected fi- 
delity is identical, J mmcm = J, lomcm , since the herald- 
ing completely removes the vacuum component. 



In general, neither measure of fidelity for either syn- 
chronized or unsynchronized systems will reach 1, except 
in the limit p — > 0. Therefore one must choose a thresh- 
old fidelity 9 that is acceptable, and then one should 
choose the largest value pe of p such that T = 0, for 
each system. Having done this, one can then compare the 
./V-fold coincidence rates. For N unsynchronized sources 
we have that p e = {2 - h - [(2 -h) 2 - 4hB} 1 ^}/2h, 
independent of N. For the same number of synchronized 
sources, p@ depends on TV and needs to be determined by 
a numerical optimisation. Figure [2] shows the resulting 
comparison of synchronized and unsynchronized systems. 
The waiting times scale exponentially with the number 
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FIG. 2: Multiphoton waiting times. The blue bars show 
the average waiting time between iV-photon events for a sys- 
tem of iV unsynchronized downconversion sources, assuming 
a pulse repetition rate of 1Z — 1 GHz, a heralding efficiency 
of h = 50% and a threshold fidelity 9 = 90%. The red 
bars show the corresponding waiting times when the system 
is synchronized with N — 1 memories, with memory efficien- 
cies r/s = T)i — 75% and a time-bandwidth product B = 1000, 
where postselection on at least iV photons is used (we set pe 
so that T = 0). The green bars show the waiting times with- 
out postselection (we set pe so that T = 0), where to achieve 
the required unpostselected fidelity threshold we now assume 
memory efficiencies of r/ s = T/ r = 99%. 



N of photons required, and without synchronization a 
12 photon experiment would require more than 30 years 
in between coincidence events, so that quantum comput- 
ing with photons using such a system is totally unfea- 
sible. However the use of memories reduces the wait- 
ing time quite dramatically. For a postselected experi- 
ment one can use inefficient memories with r\ = 56% and 
reduce the 12-fold waiting time to ~10 ms. Quantum 
memories based on Raman scattering have already been 
demonstrated with 5 > GHz [20], B > 1000 [ZEJE2] and 
7] > 50% [53] , while highly efficient and multiplexed stor- 
age in rare-earth memories is maturing and so 
these dramatic enhancements lie well within the reach of 
current technology. Without postselection more efficient 
memories are required to achieve the fidelity threshold, 
which puts the implementation beyond current techno- 
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logical capabilities, but this improved performance fur- 
ther reduces the waiting time to <1 /is. 

In summary, we have analyzed the use of quantum 
memories for the synchronization of multiple single pho- 
ton sources as a canonical application of quantum stor- 
age for the enhancement of photonic information process- 
ing. We derived an analytic formula for the multiphoton 
rate achievable and showed that the most important fig- 
ure of merit for quantum memories is the product rjB 
of the memory efficiency with its time-bandwidth prod- 
uct. Finally we extended our model to include higher- 
order photon number contributions, so that the quality 
of the states produced with and without memories could 
be compared. We showed that even inefficient memories 
can produce enormous improvements in the multiphoton 
rate when combined with postselection. Without posts- 
election, highly efficient memories are required to match 
the quality of unsynchronized sources, but if these are 
available the gain in multiphoton rate becomes larger 
still. It would be interesting to consider the effects of 
noise in the memories, or extensions to more complicated 
synchronization protocols. It is expected that similar ad- 
vantages could pertain to the scaling of other heralded 
quantum operations, such as entanglement generation or 
two-photon gates. While much attention in the quantum 
memory community has focussed on the need for long 
storage times and high efficiencies in the context of quan- 
tum repeaters [3D] , our analysis underlines the value 
of developing quantum memories for local synchroniza- 
tion, for which lower efficiencies still provide considerable 
advantages, and for which the time-bandwidth product B 
is much more important than the absolute storage time. 
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Appendix 

The readout probability R can be computed by track- 
ing the belief state y< m ) = [F (m) , V (m ^] T , where is 
the probability that we believe the memory to be charged 
at the m th time step. We have 

y (m) =Sy (m-l } . S=(® j), 

where w = Rq (z = qR) is the probability that we believe 
an empty (charged) memory becomes charged (empty) 
over the course of one time step. In the steady state 
y( m ) _ ^ y = w/(w + z). On the other hand, readout 
occurs when we believe N — 1 other photons to be avail- 
able, so we can write R — qY N ~ 2 , where Y = q + qV is 
the probability that we believe a source has provided a 
photon, either directly or through its memory. Combin- 
ing these relations we obtain the consistency condition 
(1 - 2q)Y N ~ 1 + q 2 Y N ~ 2 + Y - 1 = 0, the positive real 
root of which can be found numerically, which then fixes 
R. Note that when N = 2, we have R = q. 
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